Dynamical response function of a compressed lithium monolayer. 
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Since recent both theoretical and experimental results have proved that the simple behaviour light 
alkaline metals present at equilibrium breaks when high pressures are applied, they have become 
an important object of study in Condensed Matter Physics. On the other hand, development 
of new techniques in the atomic manipulation allows the growth of atomic monolayers (ML's), 
therefore rising the interest to analyze low dimensional systems under different conditions. In 
particular, new ab initio calculations performed for a lithium ML show that its electronic properties 
experience important modifications under pressure, which could lead to significant modifications in 
its dynamical response function. In this article we perform ab initio calculations of the dynamical 
response function of a lithium ML analyzing its evolution with increasing applied pressure. We show 
that besides the well known intraband and interband plasmons, rising electronic density induces 
characteristic features of acoustic plasmons related to the presence of two types of carriers at the 
Fermi level. 
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I. INTRODUCTION 

Light alkali metals have lately become an important ob- 
ject of study due to the complexity induced by pressure 
in their behaviour, in contrast to the simple picture they 
present under normal conditions of pressure and tem- 
perature. At equilibrium, alkalies are frequently consid- 
ered to be simple metals due to their crystallization in 
high symmetric structures, monovalency and high con- 
ductivity. At the same time the interaction between va- 
lence electrons and ionic cores, or pseudopotential, has 
been argued to be weak in such systems. As a conse- 
quence, the nearly free electron (NFE) model accurately 
describes their electronic properties in absence of applied 
pressure Q as it can be easily checked by their almost 
spherical Fermi surfaces. However, when high pressures 
are applied, the magnitude of the pseudopotential rises, 
the core electrons substantially overlap and the otherwise 
valid NFE model breaks. Neaton et al theoretically 
analyzed compressed bulk lithium, predicting that high 
pressures could induce phase transitions to less symmet- 
ric, lower coordinated structures, associated to electronic 
localizations. These theoretical predictions were exper- 
imentally confirmed by Hanfiand et al who found 
that lithium undergoes several phase transitions from a 
simple, high symmetric, bcc structure at equilibrium to 
a complex c/16 with 16 atoms per unit cell at around 
40 GPa. It is important to mention that pressure in- 
duced transitions from simple to more complex struc- 
tures are not singular to lithium but have also been ob- 
served in heavier alkalines Q. Another important fea- 



ture which shows up the complexity induced by pressure 
in the behaviour of light alkalies is that despite the su- 
perconducting transition for lithium at equilibrium has 
not been found yet, T c < 100 /^K when compressed to 
around 30 GPa, T c rises up to 15 K 0,0], becoming the 
highest transition temperature between simple elements 
U Therefore, the characterization and understanding of 
the physical properties of compressed light elements be- 
comes a priority. On the other hand, the development of 
new techniques for the atomic manipulation allows the 
growth of atomic monolayers (ML's) on inert substra- 
tum, semiconductors or noble gases. These new possi- 
bilities rise the interest to analyze physical properties of 
low-dimensional systems under different conditions, as it 
could be the case of a lithium ML under pressure. In ad- 
dition, extending our conclusions to the bulk will also give 
another perspective to understand the physical origin un- 
der the experimentally observed features in compressed 
lithium. Previous ab initio analysis of the structural and 
electronic properties of a lithium ML reveal impor- 

tant modifications in both its band structure and Fermi 
line, which will also lead to significant modifications on 
its dynamical response function, up to now just studied at 
equilibrium |ll| . In this article we perform calculations 
of the dynamical response function of a lithium ML at 
different pressures, analyzing its evolution with increas- 
ing electronic density. Besides the common intraband 
and interband plasmons, at a certain value of the lateral 
pressure applied to the ML, we observe characteristic fea- 
tures of acoustic plasmons related to the presence, at this 
pressure range, of two types of carriers at the Fermi level. 
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In Section [H] we describe the theoretical and computa- 
tional background of this work. Results and Discussion 
are presented at Section ITTll Unless otherwise stated we 
use atomic units throughout, i.e., e 2 = h = m e = 1. 



II. THEORETICAL AND COMPUTATIONAL 
BACKGROUND. 

We consider a perturbing positive charge located far from 
the ML, zq >> d (d being the thickness of the ML). 
Then, the differential cross section for a process in which 
the electron is scattered with energy lo and momentum 
transfer |q| is proportional to Im g(q, oj), where g(q;cj) 
is defined as the density response function of the ML 

pin, 

g(q;uj) = -*Lj dzj dz'x 
*XG=o,G>=o{<i.,z,z';u)e\^+*'\ (1) 

where q belongs to the two-dimensional Brillouin zone 
(2DBZ) and G and G' are reciprocal two-dimensional 
lattice vectors. In order to obtain q{c\\ui) we strictly fol- 
low the method presented in Refs. 0,0]- Thus, our 
calculations are performed within the framework of time 
dependent density functional theory (TDDFT), where 
the random phase approximation (RPA) 0| has been 
applied. We use a local density approximation (LDA) 
in order to calculate the total energy and one-particle 
eigenvalues and eigenvectors of the ground state, which 
were evaluated with the use of Troullier-Martins Li pseu- 
dopotentials 0]. The description of the ML is imple- 
mented by a supercell which contains 20 layers of vacuum 
between ML's in order to minimize interaction between 
them. Effects of compression were simulated by reducing 
the lattice parameter. For the evaluation of \%, c P re ~ 
vious step in the calculation of xg,c, we implement up 
to ^80 bands. In addition, we use ~ 7700 k-points for 
the sampling of the 2DBZ. 




FIG. 1: Band structure and density of states (DOS) for (a) a 
hexagonal lithium ML at equilibrium (r a — 3.02 a.u.) and (b) 
a square structure at r a = 2.15 a.u. (close to the structural 
transition). The Fermi energy is represented by the dotted 
line. The bandwidth of the first band decreases with pressure 
and induces the occupation of the second band. 



III. RESULTS AND DISCUSSION 

According to previous ab initio calculations, the lithium 
ML adopts an hexagonal (hex) structure at equilibrium, 
which corresponds to r s = 3.02 j^J. However, with ap- 
plying pressure, a phase transition to the more opened 
square (sq) structure is produced at r s = 2.25. In Fig. 2] 
we plot the band structure and density of states (DOS) 
of a lithium ML corresponding to the hex structure at 
equilibrium and the sq structure with 7' s = 2.15, close to 
the structural phase transition. We can observe a quasi 
free electron-like behaviour at equilibrium, where only 
one band with a parabolic-type dispersion of s-character 
is occupied. In addition, the DOS for occupied states 
is almost constant, reflecting the quasi two-dimensional 
character of the ML. However, we can notice important 
modifications at higher pressures: the first band flattens 



associated to the increasing of the band gap at the zone 
boundary, indicating a clear electronic localization. At 
the same time, as a consequence of the decrease in the 
energy difference between the first two bands at T, the 
p z -band, which is antisymmetric in the direction perpen- 
dicular to the plane of the ML, starts to be occupied. 
Therefore, we have two types of carriers at the Fermi en- 
ergy, a fact which plays a fundamental role in the pressure 
induced properties of the electronic collective excitations, 
as will be explained bellow. It is noteworthy that the 
Fermi line of the ML also suffers significant deviations 
with increasing electronic density |l0j . as it is shown in 

Fig. m 

These important pressure induced modifications in both 
the band structure and the Fermi line of the ML are 
expected to have a profound effect in the dynamical re- 
sponse function. In what follows, we will analyze the loss 
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FIG. 2: Fermi lines of a hexagonal lithium ML at equilib- 
rium, r s = 3.02 (left), and a square lithium ML at r s — 2.15 
(right). Solid line corresponds to the Fermi line of the s-band. 
For values of the density r s < 2.4, the p z -band starts to be 
occupied, so that this band is represented by a dashed line. 
The zone boundary is plotted with a dotted line. 



traband collective excitations inside the s-band, whereas 
the second one with energy of ~5 eV is related to the in- 
tcrband plasmon, which physically corresponds to collec- 
tive motion of electrons in the direction perpendicular to 
the ML. Dispersion curves for both intraband and inter- 
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function, proportional to the imaginary part of <?(q; u>), at 
two different selected densities, r s = 3.02 (equilibrium) 
and r s = 2.4. At the latter density the phase transi- 
tion has not yet been produced whereas the p z -band is 
already occupied. In such a way, we can conclude that 
the physical origin of the interesting phenomena which 
will be described bellow is not related to the hex — > sq 
phase transition, but is a direct consequence of the pres- 
ence of two different conductors at the Fermi energy. In 
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FIG. 3: Loss function Im g(q;ui) of the hexagonal lithium 
ML at equilibrium for two different values of the wavefunction 
in the TK direction. One can see one-particle electron-hole 
excitations and two peaks associated to the intraband and 
interband plasmons. 

Fig. |3 we present the loss function of a lithium ML at 
equilibrium for two selected wavevectors in the TK di- 
rection. As we are mainly interested in collective exci- 
tations, small momentum transfers are considered. Our 
spectra consists in one-particle electron-hole excitations 
with features corresponding to intraband excitations at 
small energies and interband excitations at higher ener- 
gies, together with two plasmon peaks. The first peak is 
situated at an energy close to ^3 eV and linked to the in- 
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FIG. 4: Dispersion curves of the ab initio collective intra- 
band (squares) and interband (circles) excitations in the TK 
direction. The non-zero minimum energy in the limit q — > 
indicates that electrons located at neighbouring layers feel the 
interaction. 

band collective excitations are plotted in Fig. 0] As our 
selected wavefunctions are much smaller than the char- 
acteristic Fermi momentum of the ML, we would expect 
a y/q dispersion |17| for the two-dimensional intraband 
plasmon, instead of the non-zero minimum energy for 
3D plasmons However, the Fourier transform of the 
electronic potential of electrons located at neighboring 
layers decays exponentially with |q|£, so that when deal- 
ing with vectors whose size is comparable to the inverse 
of the intermonolayer distance, q ~ 1/L = 0.05 a.u, elec- 
trons feel the interaction with other layers. Therefore, 
the plasmons show a three dimensional character at very 
small momentum, which is reflected in the gap at |q| — > 
that can be seen in Fig. 

Loss function of the ML with r s = 2.4 at selected mo- 
menta in the TK direction is presented in Fig. The 
s-bandwidth at this density is very similar to that at 
equilibrium. As a consequence, the collective excitations 
inside the s-band require a very similar energy at both 
densities, equilibrium and r s = 2.4. Thus, the intraband 
plasmon peak is situated at ~3.5 eV, close to the location 
at equilibrium. However, the energy gap between the first 
two bands, s and p z , increases under pressure, therefore 
requiring more energy the collective motions of electrons 
between these two bands to be excitated, and shifting the 
corresponding interband plasmon peak to around ^6.5 
eV. As well as the two plasmon peaks appearing also in 
absence of applied pressure, here one can also observe 
peaks associated with two other interband collective ex- 
citations with energies ~7.7 eV and ~9.5 eV, due to the 
presence of two partly occupied energy bands. In addi- 
tion, a remarkable peak is induced by pressure at energies 
close to zero. Besides the common intraband and inter- 
band collective excitations, the presence of two type of 
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FIG. 5: Loss function Im g(q; uj) of the hexagonal lithium ML 
with r s = 2.4 for two different values of the wavefunction in 
the TK direction. Besides the common intraband plasmon, 
close to the same energy as at equilibrium, and interband 
one, located at a higher energy, one can also notice a remark- 
able peak at small uj, the characteristic feature of an acoustic 
plasmon. 

carriers at the Fermi level with different velocities allows 
the existence of another kind of collective modes, known 
as acoustic plasmons because of their linear dispersion in 
|q| 0|, as it can be checked in Fig. HJ1 The same acoustic 



plasmon was also predicted for a berilium ML [llj . Con- 
trary to the case of two dimensional- like plasmons, which 
are well established experimentally jl^, 120] , the acoustic 
counterparts have not yet been unambiguously identi- 
fied. It is known that, in this case, the light electrons 
screen the Coulomb repulsion between the electrons with 
a lower velocity, so that an additional intraband acoustic 
collective excitation is produced. 
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FIG. 6: Dispersion curve of the ab initio acoustic plasmon 
induced by the applied pressure in a hexagonal lithium ML 
with r s = 2.4. A clear lineal dependence with |q| can be 
observed. Dots are calculations whereas dashed line is a fit. 
The slope of the fit is close to the velocity of the slow electrons. 
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